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Abstract 

The nondetection of neutrinos coming from Gamma Ray Bursts (GRBs) by the IceCube exper- 
iment has raised serious questions on our understanding of GRB's and the mechanism of neutrino 
flux production in them. Motivated by this and the need for a precise calculation for GRB neutrino 
flux, here we study the effects of beyond standard model physics on the GRB neutrino flux. In 
the internal shock model of GRB, high energy neutrinos are expected from muon, pion and kaon 
decays. Using the latest best fit neutrino oscillation parameters, we compute the expected fiux on 
earth for standard as well as non-standard oscillation scenarios. Among the non-standard scenar- 
ios, we consider neutrino decay, pseudo-dirac nature of neutrinos and presence of one eV scale light 
sterile neutrino. Incorporating other experimental bounds on these new physics scenarios, we show 
that neutrino decay scenario can significantly alter the neutrino fiux on earth from the expected 
ones whereas the corresponding changes for pseudo-dirac and sterile neutrino cases are moderate. 
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I. INTRODUCTION 



Recent neutrino oscillation experiments have provided significant amount of evidence 
which confirms the existence of the non-zero yet tiny neutrino masses [1-5]. The smallness 
of neutrino masses compared to electroweak scale can naturally be explained by seesaw 
mechanism [6-10], the simplest version of which corresponds to the inclusion of three singlet 
right handed neutrinos into the standard model. Although these seesaw models can naturally 
explain the smallness of neutrino mass compared to the electroweak scale, we still do not 
have a complete understanding of the origin of neutrino mass hierarchies as suggested by 
experiments. Recent neutrino oscillation experiments T2K [11], Double ChooZ [12], Daya- 
Bay [13] and RENO [14] have not only made the earlier predictions for neutrino parameters 
more precise, but also predicted non-zero value of the reactor mixing angle 6*13. The latest 
global fit value for 3cr range of neutrino oscillation parameters [15] are as follows: 

Am^^ = (7.00 - 8.09) x 10'^ eV^ 

Aml^ (NH) = (2.27 - 2.69) x 10"^ eV^ 
Aml^ (IH) = (2.24 - 2.65) x 10"^ eV^ 

sin^^ia = 0.27-0.34 

sin^^ss = 0.34 - 0.67 
sin^^ia = 0.016 - 0.030 (1) 

where NH and IH refers to normal and inverted hierarchy respectively. 

Although the next generation neutrino oscillation experiments are expected to shed light 
on the origin of mass hierarchies as well as the Dirac CP phase, it is worth exploring if there 
exists an alternate experimentally verifiable way to understand some of the yet unresolved 
issues in neutrino physics. It will be even more exciting if such alternate ways can also 
confirm or rule out some of the well motivated beyond standard model frameworks which 
may or may not be seen in collider experiments. It turns out that the neutrino telescopes 
which have been designed to observe high energy cosmic rays, can be a promising setup to 
search for new physics. 

The consequences of many such well motivated new physics scenarios on the observations 
of neutrino flux observed by neutrino telescopes have been studied by several groups [16-31]. 
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Motivated by these, here we pursue a similar study on the possibihty of observing new physics 
at neutrino telescopes. In particular, we focus on high energy neutrinos coming from Gamma 
Ray Bursts (GRBs) and present an analysis of how the expected total flux of neutrinos at 
neutrino telescopes can change significantly by the presence of new physics. Among new 
physics scenarios, we consider neutrino decay, presence of one light sterile neutrino and 
pseudo dirac nature of neutrinos. 

Inside a ORE, shock accelerated protons may interact with low energy photons leading 
to the production of high energy mesons p7 — ?■ vr+'^X and these pions subsequently decay 
to high energy neutrinos Tr"*" — /i'^i^^, fi~^ — )■ e'^u^Uf.. This flux has already calculated in 
[32-42]. Recently IceCube collaboration has claimed to reach the sensitivity of detecting 
neutrino flux from GRBs at TeV energy [43]. The combined operation of Icecube 40 and 
59 string for the time period of April 5, 2008 to May 2010 for GRB neutrinos has placed 
a tighter upper bound, 3.7 times below the theoretical predictions [44]. [45-47] have re- 
calculated the neutrino flux from the 215 GRBs used by IceCube during their period of 
detection, and concluded that the neutrino flux predicted theoretically in the papers by 
IceCube collaborations is an overestimation. 

However p-f in GRBs can also produce high energy neutrons, and they will decay as 
n p + e~ + to antineutrinos [48]. The other secondary products in py interactions 
are p7 — K~^'^X where X can be either A°, S'^ or Kaons decay to lighter mesons, 
leptons and neutrinos [49-51]. In [49, 50] the dominant decay channel kaon to neutrino, 
— > 1/^(63%) was taken while we have considered all the channels of K~^'^ decaying to 
neutrinos [5]. So one can find the total neutrino flux from GRBs has a contribution from 
these processes too. 

Using the above mentioned possible origin of high energy neutrino flux from GRBs and the 
best fit values of neutrino oscillation parameters (1), we show that the individual neutrino 
flux can change significantly from the ones expected from standard oscillation paradigm. 
More specifically, the scenario of neutrino decay can change the expected flux to a great 
extent, wheareas the changes in the scenario of sterile neutrino and pseudo-dirac neutrino 
are somewhat moderate and should be detectable in future neutrino telescopes. 

This paper is organized as follows: In section II, we discuss the possible ways high energy 
neutrinos can originate in GRBs. Then we discuss the change in neutrino flux due to 
standard as well as non-standard oscillations in section III and finally conclude in section 
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IV. 



II. HIGH ENERGY NEUTRINOS FROM GAMMA RAY BURSTS 

In calculating the high energy neutrino flux from individual GRBs we have used the 
method as in ref. [51]. Frames of references are assigned as "c" for comoving or wind rest 
frame, "p" for proton rest frame. Quantities measured in the source rest frame are written 
without any subscript. Shock accelerated high energy protons interacting with low enrgy 
photons will produce high energetic muons, pions, neutrons and kaons. We have used the 
low energy photon flux typically observed by Swift in the energy range of 1 KeV to 10 
MeV to calculate the neutrino flux from individual GRBs with break at in the source rest 
frame related to the break energy in the comoving frame ^ as e!^ = Fe!^ ^. 

7i < 2, and 72 > 2. The normalization constant A is related to the internal energy density 
Uhj, 



71-2 



L72-2 71-2J 

The maximum energy of the shock accelerated protons in the GRB fireball can be calculated 
by comparing the minimum of the py interaction time scale (tp^), p-synchrotron cooling time 
scale (tsyn) and dynamical time scale (tdyn) of a GRB with the acceleration time scale (tacc) 
of the protons as discussed in [52]. 

^acc Tflinitpy, tdyni ^syn) (3) 

We have considered production of pions in p-f interactions with two particle final states 
through the decay of resonant particle A"*". At the delta resonance both 7r° and vr"*" have 
been assumed to be produced with equal probabilities, vr"^ gets on the average 20% of the 
proton's energy. The charged pions decay to muons and neutrinos. Finally the muons decay 
to electrons and neutrinos, antineutrinos. Each pion decay followed by muon decay gives two 
neutrinos, one antineutrino and one positron. If the final state leptons share the pion energy 
equally then each neutrino carries 5% of the initial proton's energy. The specific parameters 
of GRBs can be denoted as the following way, the fireball Lorentz factor F300 = F/300, 



photon luminosity ^^7,51 = L^^j {IQi^^ergs /sec), variability time t^^^z = {t^llQ~'^sec) are the 
important parameters of a GRB. The internal energy density U relates to photon luminosity, 
= Attt^'T'^cU . = r'^cty is the internal shock radius. 

The total energy to be emitted by neutrinos of energy e,^ from photo-pion decay (con- 
sidering muon and pion decay neutrinos together) in the source rest frame of a GRB is 
[38], 



, dK^ 3/. (1 - 6e - es) I 1 e,,. < et^^ 



where El^^° is the total isotropic energy of the emitted gamma-ray photons in the energy 
range of IkeV to lOMeV, which is available from observations. It is the product of with 
the duration of the prompt emission from the GRB. and es are the energy fractions carried 
by electrons and the magnetic field respectively. The maximum energy of neutrinos from 
pion decay is approximately 5% of the maximum energy of protons (ep^max)- ^t,ij. presents 
the neutrino energy where muon synchrotron cooling starts. 

Ultrahigh energy neutrons are also produced in p7 interactions along with pions and 
kaons. These neutrons (with Lorentz factor r„) decay (n — )■ p + e~ + z7e) to with a decay 
mean free path cTnfn = 10{en/EeV) Kpc. r„ = 886 seconds is the lifetime of a neutron 
in its rest frame and en is its energy in the source rest frame. Assuming the probability of 
production of neutrons in resonant p7 interactions to be half one would be able to calculate 
the fraction of a proton's energy lost to neutron production in the process p7 — ?■ ir^n at the 
A— resonance [48] as, 

7i"l ^ 



where fn = 4 ^'^^^'f . 1 1- and ^„ = 0.8. And the energy flux of antineutrinos 

-I- 300 *i'>-3'^7,A/eV \ —^r ^\ 

L72-2 71-2J 

of energy ep,„ can be estimated 
travelling a distance Dg [53] as. 



L72-2 71-2J 

of energy ep,„ can be estimated with the neutron flux {dNn/ den) at the source rest frame 



2 dNp 



den dNn Dsrnn \ m„ 



1 -e 



en den ^ ^ 2 eo 



^0 



X (6) 



eo is the mean energy of an antineutrino in the neutron rest frame. 
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In p7 interactions charged kaons {K~^) are produced through the following interactions 
P7 — )■ K~^A^ and p-f — )■ K~^T,^. Although the cross-sections of these interactions are lower 
compared to photo-pion production [54], at very high energy the total neutrino flux produced 
through kaon decay becomes higher than that from photo-pion decay [49, 55]. The fractional 
energy transferred from shock accelerated protons to kaons has been calculated for all the 
channels of from both resonant and multiparticle production from p7 interaction as [51]. 

decay to secondary neutrinos by the following channels, — )■ yU"*"i^^(63%), 7r+7r*^(21%), 
7r+7r+7r-(6%), 7r°e+z/e(5%), 7i'^fi+Uf,{3%) and 7i+7i°n%2%). Due to their heavier mass K+ 
cools at higher energy compared to muon and pion. The synchrotron cooling break energy 
in the kaon spectrum is at 

6^ = 2.2 X lO^e^^e^^/^L^Jf GeV. (7) 

It is derived by comparing the decay and cooling time scales of kaons. Although the cross- 
section of kaon production is much less than that of pion production, the neutrino flux from 
kaon channel exceeds the flux from pion channel at very high energy due to the slower rate 
of cooling of kaons. 

Neutral kaons are produced in py interactions with a cross-section ^ ^ 0.6 x 10^^^ 
cm^ [56] at the peak energy eo^" = 1.45 GeV and width 5e/<o = 0.7 GeV. Half of the neutral 
kaons are assumed to be long lived kaons (-R'^). can be produced in interactions 
with multiparticle final states {jyy K'^A^ir^, i^'£A°7r+ and K^T.+ir'^). The cross-sections 
of these interactions are measured as 0.5 x 10~'^°cm^, 0.5 x lO^^'^cm^ and 0.2 x lO^'^^cm^ 
respectively [56]. decays through the following channels 7r~'^e~z7g(39%), 7r+/i^z7^(27%), 
^o^o^o^21%)^ and 7r"''7r~7r°(13%). Kg decays to two charged pions through this channel 
7r+7r~(69%). The pions finally decay to neutrinos and antineutrinos. The total neutrino 
flux from and Kg can be calculated in the same way as K~^ . 

The observed total neutrino flux with energy on earth is, 

where z is the redshift of the GRB. 

Our calculations are based on the standard internal shock model of GRBs. In internal 
shocks the shock radius is related to the bulk Lorentz factor F and variability time t^, 
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Td = r^ct^, where c is the speed of hght. We have not assumed any relation among the GRB 
parameters F and isotropic energy [57-59] or peak luminosity and observed break energy in 
the low energy photon spectrum [59]. We have taken the following 4 set of GRB parameters 
to calculate the neutrino flux. 

1. 7i = 1, 72 = 2.2, = 10^3 erg/sec, T = 600, ty = 20 msec, = 0.5MeV, e^/ = 1, 
/° = 0.09 and = 2.16 x 10^^ cm. 

2. 71 = 1.2, 72 = 2.5, = 10^^ erg/sec, T = 600, = 20 msec, = 0.5MeF, e^/ 
ee = 10, {ee = 0.6, = 0.06), /° = 0.17 and = 2.16 x 10^^ cm. 

3. 7i = 1.8, 72 = 2.01, = 5 X 10^^ erg/sec, T = 130, t„ = 25 msec, = 0.5Mel^, e^/ 

= 1, /o = 0.16 and = 1.26 x lO^^ ^ni. 

4. 7i = 1.2, 72 = 2.2, = 10^^ erg/sec, T = 1000, = 20 msec, = 0.5Me1/, e^/ 

= 1, /o = 0.12 and = 6 x 10^^ cm. 

III. EFFECT OF NEW PHYSICS ON GRB NEUTRINOS 
A. Standard Neutrino Oscillation 

Neutrino oscillation data clearly indicate the smallness of three Standard Model neutrino 
masses [1-4] which can be naturally explained via see-saw mechanism [6-9]. Without using 
any particular type of seesaw, here we use the most general neutrino mixing MNS matrix 
[60] and the standard vacuum oscillation probability given by 

P(z/« ^ Up; L) = 6^p-J2 U:^Up,U^,U;,il - e^^^^-^) (9) 

where Uai is an element of the the MNS matrix, a, i denoting flavor and mass eigenstates 
respectively. For cosmological distances like the typical distance of a GRB from earth we 
can assume the limit L — t- oo which simplifies the above expression for probability to 

PiUa^Up) = J2\Uaj\'\Up,\' (10) 
j 
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Yj=1.2, 72=2.5, L^IO erg/sec, r=600, t^=20 msec, 
t|,''=0.5 MeV, eg/ 6^=10. 
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FIG. 1: Total muon neutrino flux on earth for standard oscillation i.e., without decay (labelled as 
WOD) and for neutrino decay scenario in NH regime 

The role of standard neutrino oscillation on ultra high energy neutrino flux from GRB's were 
studied earlier in [61, 62]. These studies concluded that the neutrino flavor ratio 1 : 2 : at 
source would reduce to 1 : 1 : 1 at earth due to oscillations. The neutrinos while coming from 
a distance of z will undergo oscillation. We have taken the standard oscillation parameters 
along with the recent calculated value of sin^26'i3 = 0.1, where ^13 is the neutrino mixing 
angle measured by Double ChooZ [12], Daya-Bay [13] and RENO [14] collaborations and 
other oscillation parameters from global fit data [15]. 
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FIG. 2: (a) Electron, muon and tau neutrino flavor composition for standard oscillation (WOD) 
and for neutrino decay scenario (WD) for rs/ma = T2/m2 = Ti/m\ = 100 for GRB 1. (6) Neutrino 
flavor flux ratio for the GRB 1 

B. Neutrino Decay 

If the neutrino mass eigenstates are hierarchical, then a higher mass eigenstate can decay 
into a lower mass eigenstate. The role of such neutrino decay on neutrino flavor flux was 
studied in [16]. Here we consider the simplest possible situation where the heavier mass 
eigenstate completely decays into the lightest mass eigenstate which is kinematically stable. 
Thus in case of normal hierarchy (NH)mi,3 > rriy^ > m^-^, the mass eigenstate ratio at earth 
will be 1 : : 0. Hence the flavor ratio at earth will be 0.67 : 0.26 : 0.07 for NH . 

For incomplete decay, we have to include the decay factor in the expression for probability 
of oscillation. This decay factor which accounts for the depletion in neutrino flux due to 
the decay of mass eigenstate rrii with rest-frame lifetime Tj and energy E propagating over 
a distance L, is exp(— ^^). The expression for neutrino flux at earth in this case becomes 

<p,^{E) = 5^5](/)j(^)|[/;,,n[/„,iv-- (11) 

where 0° denotes the flux of neutrino flavor /3 at source. Taking the flavor ratio at source 
to he X : y : z and using the oscillation parameters as in [15], we compute the flux for 
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FIG. 3: Total muon neutrino flux on earth for standard oscillation (SO) and for Pseudo dirac 
neutrino scenario 

different neutrino flavors at earth. For the purpose of our calculation we consider three 
cases Ts/ms = Tilra2 = Ti/nii = 0.01, 10, 10^ s/eV respectively as shown in figure 1. Also, 
the electron, muon and tau flavor composition of neutrino flux on earth is shown for one 
speciflc choice of GRB parameters as well as r/m in flgure 2 (a). We also show the flavor 
ratio of muon type to electrom plus tau type neutrinos in flgure 2 (b) for the same choice of 
GRB parameters. 
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FIG. 4: (a) Electron, muon and tau flavor composition of neutrino flux on earth for standard 
osciflation (SO) and for Pseudo dirac neutrino scenario (PD) 2 — 10"^^ for GRB 1. (5) 
Neutrino flavor flux ratio for the GRB 1 

C. Pseudo Dirac Neutrinos 

Instead of a being purely Dirac with a mass m^), neutrino can be a mixture of two almost 
degenerate Majorana neutrinos. Please see [21] for the role of pseudo Dirac neutrinos in ultra 
high energy neutrino flux and references therein for earlier works on pseudo Dirac neutrinos. 
In such a case the majorana mass terms mj^, m/j <^ m/) and the the mass splitting is ~ 
2mo{mL+'m{i). Thus the neutrino mass and flavor basis become (z/^, 1/2 ,1^^ ,1^1 ,1^2 ^'^3) 
(z/g, z/^, z/t-, i/g, u'^, u'^) respectively. As shown in [21], the neutrino flavor conversion probability 
can have a form as simple as 

Po^p = (12) 

a j=l 

The new constribution coming from pseudo Dirac nature of neutrinos will be negligible until 
E/L becomes of the order of Sm'j. 5m? can be as large as 10"^^eV^ for z/1^2 and as large 
as 10~^eV^ for z/3. Taking the initial neutrino flavor ratio at source to he x : y : z and 
using the same neutrino oscillation parameters as in [15] we find the probability of detecting 
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FIG. 5: Total muon neutrino flux on earth for normal oscillation (SO) and for one extra sterile 
neutrino (SN) case 

individual flavors on earth. Here we use 5m1 2 = 10"^^, 10~^^ eV^, 5m| = 10~^ eV^ for the 
purpose of our calculation. L can be taken to be lOOMpc which is the typical distance of a 
GRB from earth. The consequences of this pseudo-dirac nature of neutrinos on the muon 
neutrino flux on earth is shown in figure 3. The electron, muon and tau flavor composition 
of neutrino flux for one specific GRB parameter and is shown in figure 4 (a). We also 
show the flavor ratio of muon type to electrom plus tau type neutrinos in figure 4 [b). 
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FIG. 6: (a) Electron and muon flavor composition of neutrino flux on earth for normal oscillation 
(SO) and for one extra sterile neutrino (SN) case for Si = 0.5 for GRB 1. (6) Neutrino flavor flux 
ratio for the GRB 1 

D. Presence of Sterile Neutrinos 

Precision measurement of the Z boson decay width restricts the number of standard model 
neutrinos to three. However we can still have an eV scale neutrino which has no coupling 
to the Z boson and hence called sterile. The presence of such a sterile neutrino is also 
compatible with the recent Wilkinson Microwave Anisotropy Probe (WMAP) collaboration 
data [63] with 95% allowed mass range < 0.48 eV [64]. 

The effect of sterile neutrino on ultra high energy neutrino flux was studied earlier in 
[20]. We follow their approach and use the recent oscillation data [15] to see the effect of 
such sterile neutrinos on the neutrino flux coming from GRBs. As considered in [20], we 
also consider the presence of a sterile neutrino which is almost degenerate with the lightest 
active neutrino mass eigenstate. Interestingly, if the mass difference is (5m^ < 10~^^ GeV 
then there is no experimental constraint on the mixing angles of this sterile neutrino with 
the three active ones. As outlined in [20], the simplified oscillation probability for 3 + 1 
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neutrino scheme is 

-Pa/3 = 'S'i|t/Q,ip|f//3i|^ + S2\Ua2\^\Ui32\^ + 5*3 1 f/^s ^ | t/gs | ^ (13) 

where Si = cos^ (f>i + sin^ 0j and (pi is the mixing angle between the active neutrino mass 
eigenstate Ui and the sterile state Us. In general Si factors can vary in the range | < 5^ < 1. 
Taking the flavor ratio at source to he x : y : z and using the neutrino oscillation data 
from [15] we find the probability of individual neutrino flavors at earth. For zero mixings 
that is, (pi = 0, the flavor ratios correspond to the standard oscillation as discussed earlier. 
To show the extent to which the presence of sterile neutrino can affect the flavor fluxes we 
plot the total muon neutrino flux on earth for two different cases Si = 0.5, 0.8 and compare 
with the standard predictions as shown in figure 5. The maximum value of this parameter 
S = 1 reproduces the standard oscillation results as seen from figure 5. We also show the 
electron, muon and tau flavor composition of neutrino flux on earth in figure 6 (a). We 
also show the flavor ratio of muon type to electrom plus tau type neutrinos in figure 6 
(b). It should be noted that a larger (~ eV^) mass difference between the active and the 
sterile mass eigenstate has other interesting motivations like explaining the anomalies found 
in MiniBooNE [65] and LSND [66] which we do not attempt to explain here. In case of 
larger mass difference, the mixing angles between active and sterile neutrinos will be tightly 
constrained and hence will not affect the flavor ratios of ultra high energy neutrinos. 

IV. RESULTS AND DISCUSSION 

We have presented an analysis of the effect of three different new physics scenarios namely, 
neutrino decay, pseudo Dirac nature and presence of one light sterile neutrino on the flux of 
high energy neutrinos coming from GRB's. We first calculate the GRB neutrino flux from 
various sources like pion, muon, neutron and kaon decays by choosing four sets of GRB 
parameters. We then incorporate the standard neutrino oscillation between all three flavors 
(electron, muon and tau) and calculate the total muon flux on earth as well as the flavor 
composition. We show that different choices of GRB parameters can give rise to different 
neutrino flux on earth. 

After calculating the neutrino flux on earth using standard neutrino oscillation physics, 
we then incorporate the three beyond standard model frameworks mentioned above to calcu- 
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late the total muon flux on earth as well as the flavor composition. We flnd that for neutrino 
decay scenario, the changes in neutrino flux from the standard oscillation case can be very 
signiflcant and could give rise to an explanation of the present non-detection of GRB neutri- 
nos at IceCube experinet. However for other two beyond standard model physics scenarios 
we consider namely, for pseudo Dirac and sterile neutrino cases, the changes are moderate. 

For the standard oscillation case, we get three plateau regions in the flavor ratio plot as 
can be seen in part (b) of flgure 2, 4 and 6. At lower energy, the neutron decay channel 
dominates while the muon damping corresponds to the transition to the second plateau 
region. As can be seen from the same flgures, for the three beyond standard model scenarios 
the flavor ratio plot changes from the standard oscillation case. 

V. ACKNOWLEDGEMENT 

We would like to thank Prof. Nayantara Gupta, RRI Bangalore and Prof. S. Umasankar, 
IIT Bombay for comments and suggestions. 



S. Fukuda et al. (Super-Kamiokande), Phys. Rev. Lett. 86, 5656 (2001), hep-ex/0103033. 

Q. R. Ahmad et al. (SNO), Phys. Rev. Lett. 89, 011301 (2002), nucl-ex/0204008. 

Q. R. Ahmad et al. (SNO), Phys. Rev. Lett. 89, 011302 (2002), nucl-ex/0204009. 

J. N. Bahcall and C. Pena-Garay, New J. Phys. 6, 63 (2004), hep-ph/0404061. 

K. Nakamura et al. (Particle Data Group), J.Phys. G37, 075021 (2010). 

P. Minkowski, Phys. Lett. B67, 421 (1977). 

M. Gell-Mann, P. Ramond, and R. Slansky (1980), print-80-0576 (CERN). 

T. Yanagida (1979), m Proceedmgs of the Workshop on the Baryon Number of the Universe 

and Unified Theories, Tsukuba, Japan, 13-14 Feb 1979. 

R. N. Mohapatra and G. Senjanovic, Phys. Rev. Lett. 44, 912 (1980). 

J. Schechter and J. Valle, Phys.Rev. D22, 2227 (1980). 

K. Abe et al. (T2K Collaboration), Phys.Rev.Lett. 107, 041801 (2011), 1106.2822. 

Y. Abe et al. (DOUBLE-CHOOZ Collaboration), Phys.Rev.Lett. 108, 131801 (2012), 

1112.6353. 



15 



F. An et al. (DAYA-BAY Collaboration), Phys.Rev.Lett. 108, 171803 (2012), 1203.1669. 

J. Ahn et al. (RENO collaboration), Phys.Rev.Lett. 108, 191802 (2012), 1204.0626. 

M. Gonzalez-Garcia, M. Maltoni, J. Salvado, and T. Schwetz (2012), 1209.3023. 

J. F. Beacom, N. F. Bell, D. Hooper, S. Pakvasa, and T. J. Weiler, Phys. Rev. Lett. 90, 

181301 (2003), hep-ph/0211305. 

J. F. Beacom, N. F. Bell, D. Hooper, S. Pakvasa, and T. J. Weiler, Phys.Rev. D68, 093005 
(2003), hep-pli/0307025. 

A. Bhattacharya, S. Choubey, R. Gandhi, and A. Watanabe, Phys.Lett. B690, 42 (2010), 
0910.4396. 

A. Bhattacharya, S. Choubey, R. Gandhi, and A. Watanabe, JCAP 1009, 009 (2010), 
1006.3082. 

P. Keranen, J. Maalampi, M. Myyrylainen, and J. Riittinen, Phys. Lett. B574, 162 (2003), 
hep-ph/0307041. 

J. F. Beacom et al., Phys. Rev. Lett. 92, 011101 (2004), hep-ph/0307151. 

D. Hooper, D. Morgan, and E. Winstanley, Phys.Lett. B609, 206 (2005), hep-ph/0410094. 

D. Hooper, D. Morgan, and E. Winstanley, Phys.Rev. D72, 065009 (2005), hep-ph/0506091. 

D. Meloni and T. Ohlsson, Phys.Rev. D75, 125017 (2007), hep-ph/0612279. 

Z.-z. Xing and S. Zhou, Phys.Lett. B666, 166 (2008), 0804.3512. 

A. Esmaih and Y. Farzan, Nucl.Phys. B821, 197 (2009), 0905.0259. 

R. Cowsik, T. Madziwa-Nussinov, S. Nussinov, and U. Sarkar, Phys.Rev. D86, 045024 (2012), 
1206.0713. 

D. Meloni and T. Ohlsson, Phys.Rev. D86, 067701 (2012), 1206.6886. 

P. Baerwald, M. Bustamante, and W. Winter, JCAP 1210, 020 (2012), 1208.4600. 

A. Esmaih and Y. Farzan (2012), 1208.6012. 

S. Pakvasa, A. Joshipura, and S. Mohanty (2012), 1209.5630. 

M. Vietri, Astrophys.J. 453, 883 (1995), astro-ph/9506081. 

E. Waxman, Phys. Rev. Lett. 75, 386 (1995), URL http://link.aps.org/doi/10.1103/ 
PhysRevLett . 75 . 386. 

J. Alvarez-Muniz and F. Halzen, Phys.Rev. D63, 037302 (2001), astro-ph/0007329. 
N. Gupta, Phys.Rev. D65, 113005 (2002), astro-ph/0201509. 

D. Guetta, D. Hooper, J. Alvarez-Muniz, F. Halzen, and E. Reuveni, Astropart.Phys. 20, 429 



16 



(2004) , astro-ph/0302524. 

[37] P. Bhattacharjee and N. Gupta, Phys.Rev.Lett. (2005), hep-ph/0501191. 

[38] N. Gupta and B. Zhang, Astropart.Phys. 27, 386 (2007), astro-ph/0606744. 

[39] K. Murase and S. Nagataki, Phys.Rev. D73, 063002 (2006), astro-ph/0512275. 

[40] J. K. Becker, W. Rhode, P. L. Biermann, and K. Muenich, Astropart.Phys. 28, 98 (2008), 

astro-ph/0607427. 
[41] E. Waxman, Science 315, 63 (2007), astro-ph/0701168. 

[42] K. Murase, K. loka, S. Nagataki, and T. Nakamura, Phys.Rev. D78, 023005 (2008), 0801.2861. 
[43] R. Abbasi et ah (IceCube Cohaboration), Phys.Rev.Lett. 106, 141101 (2011), 1101.1448. 
[44] R. Abbasi et ah (IceCube Collaboration), Nature 484, 351 (2012), 1204.4219. 
[45] S. Hummer, P. Baerwald, and W. Winter, Phys.Rev.Lett. 108, 231101 (2012), 1112.1076. 
[46] Z. Li, Phys.Rev. D85, 027301 (2012), 1112.2240. 

[47] H.-N. He, R.-Y. Liu, X.-Y. Wang, S. Nagataki, K. Murase, et al., Astrophys.J. 752, 29 (2012), 
1204.0857. 

[48] R. Moharana and N. Gupta, Phys.Rev. D82, 023003 (2010), 1005.0250. 

[49] P. Baerwald, S. Hummer, and W. Winter, Phys.Rev. D83, 067303 (2011), 1009.4010. 

[50] P. Baerwald, S. Hummer, and W. Winter, Astropart.Phys. 35, 508 (2012), 1107.5583. 

[51] R. Moharana and N. Gupta, Astropart.Phys. 36, 195 (2012), 1107.4483. 

[52] N. Gupta, Astropart.Phys. 31, 359 (2009), 0903.4733. 

[53] L. A. Anchordoqui, H. Goldberg, F. Halzen, and T. J. Weiler, Phys.Lett. B593, 42 (2004), 
astro-ph/0311002. 

[54] K. Glander, J. Barth, W. Braun, J. Hannappel, N. Jopen, et al, Eur.Phys.J. A19, 251 (2004), 
nucl-ex/0308025. 

[55] K. Asano and S. Nagataki, Astrophys.J. 640, L9 (2006), astro-ph/0603107. 

[56] R. Lawall, J. Barth, C. Bennhold, K.-H. Glander, S. Goers, et al, Eur.Phys.J. A24, 275 

(2005) , nucl-ex/0504014. 

[57] E.-W. Liang, S.-x. Yi, J. Zhang, H.-J. LV, B.-B. Zhang, et al., Astrophys.J. 725, 2209 (2010), 
0912.4800. 

[58] J. Lu, Y.-C. Zou, W.-H. Lei, B. Zhang, Q.-W. Wu, et al., Astrophys.J. 751, 49 (2012), 
1109.3757. 

[59] G. Ghirlanda, G. Ghisellini, R. Salvaterra, L. Nava, D. Burlon, et al. (2012), 1210.1215. 

17 



[60] Z. Maki, M. Nakagawa, and S. Sakata, Prog. Theor. Phys. 28, 870 (1962). 

[61] J. G. Learned and S. Pakvasa, Astropart. Phys. 3, 267 (1995), hep-ph/9405296. 

[62] H. Athar, M. Jezabek, and O. Yasuda, Phys. Rev. D62, 103007 (2000), hep-ph/0005104. 

[63] E. Komatsu et al. (WMAP), Astrophys. J. Suppl. 192, 18 (2011), 1001.4538. 

[64] J. Hamann, S. Hannestad, G. G. Raffelt, I. Tamborra, and Y. Y. Y. Wong, Phys. Rev. Lett. 

105, 181301 (2010), 1006.5276. 
[65] A. Aguilar-Arevalo et al. (The MiniBooNE Collaboration), Phys.Rev.Lett. 105, 181801 (2010), 

1007.1150. 

[66] A. Aguilar et al. (LSND Collaboration), Phys. Rev. D64, 112007 (2001), hep-ex/0104049. 



18 



